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INTRODUCTION 


In most practical cases the noise sources of an electric circuit are stationary 
with time, or the changes are so slow that the problem can be treated as stationary. 
There occur, however, cases where the variations in shot noise may be fairly 
fast. 

As an example, take the shot noise generated in a temperature-limited diode 
operating on an A.C, filament supply. Here the probability of electron emission 
varies periodically with a fundamental frequency equal to twice the supply line 
frequency. Even in this case, as in the stationary shot noise, the emission of an 
electron occurs independently of the number of electrons emitted before, or the 
number of those that will be emitted after the event, but the probability of emission 
is not constant with an arbitrary shift of time. 

Another example of the nonstationary shot noise can be observed in the ionization 
current of a chamber if the intensity of the radiation entering the chamber is varied 
with time, The arrival and absorption of each quantum and consequently the 
occurrence of the current pulse it produces is (when the space charge is negligible) 
independent of other quanta, but the probability of the current pulse is proportional 
to the radiation intensity and accordingly varies with time. 

Expressions for the autocorrelation function and the statistical power spectrum 
of the nonstationary shot noise are derived in this paper. Special attention is given 
to the periodic case, where the probability of current pulses is a periodic function 
of time. 
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PROBABILITY DISTRIBUTIONS 


Let us denote by h (t) the non-negative function, which determines the probability 
density of the current pulses at the moment t. Thus, the probability p (1, t) dt of 
the event that one pulse falls in the interval (t, t + dt) is 


(1) p (1, t) dt =h (t) dt 


The probability of more than one pulse falling in the infinitesimal interval is 
of a lower order of magnitude and can be considered to be zero, Hence, the only 
possible event in addition to (1, t) is (0, t) or no pulses in the interval, Thus 


(2) p (0, t)dt+p(l, t)dt=1 
From equations (1) and (2) we get 


(3) p (0, t) dt =1-h(t)dt 


In order to find the probability P (0, t,, t.) of the event that no pulses occur 
in the finite time interval (tz, to), we divide this interval into M sections, each 
of duration At =(t» - t,)/M, which is so short that h (t) does not change appreciably 
during it. Then we may write 


(4) p(0, t) At=1-h(t) At =exp[-h(t) At] 


The probability of no pulses during the time from t, to t, is equal to the joint 
probability of no pulses in any of the intervals At. These events are statistically 
independent and their joint probability is the product of the individual probabilities, 


P(,t.t)) = exp[-h(t, +m 4] 
m 
(5) 


M-1 
= exp[ -h(t; +m At)] 
m=0 


Now, as At-—>0, the approximation in equation (4) becomes exact, *) and the 


t 
2 
sum in the exponent (5) approaches the integral - J h (x) dx and 
(6) P (0, ty» to) = exp [ - [ h (x) dx | 
t 
1 


Next, to calculate the probability P (N, t;, to) of N pulses occurring in the 
interval (t), to) we consider the event as consisting of two disjoint events: 


a) In the interval (t,, tp- At) there are N-1 pulses and between At and ty 


there is one pulse, 
b) In the interval (t,,to- At) there are N pulses and between ty- At and ty no 


pulses, 
The probabilities P, and P}, of events a) and b) are evidently 


Pa (N, ty, to) P(N-1, ty, tg- At) p(1, to- At) At 


(7) 


P, (N, to) P(N, tg - At) At) At 


*) This derivation is not very rigorous, but gives a correct result. . 


| 


Since the events are disjoint, 


P(N, to) (Ne ty, tp) + Ph (N, ty» ta) 


or by equations (1) and (3) 


P(N, ty, to) = P(N-1, ty, At) - At) At 


+ P(N, ty, ty - At) An) At] 


which may be written 


P(N, ty, to) - P(N, ty, &- At 
At 


+ h(ty - At) P(N, ty, tg - At) = 


h(t, - At) P(N-1, ty, tg- At) 


and, when the limit At -0 is taken, we have the differential equation x) 


(8) P(N, ty, to) + P(N, ty, to) = h(tg) P(N-1, ty, tg) 


With the initial condition P(N, t,, tj) = 0, equation (8) has the solution 


tg tg tg 
(9) P(N, ty, = exp [- [ hexydx] fropPav-t exp [ dz] dy 
ty 


This is a recursion formula by which P(N, t,, tp) can be obtained from P (N-1, ty» ty). 
By starting from (6) N=0 one can easily arrive at the result 


*) The derivation of (8) is analogous to that given by W.B.Davenport Jr. and 
W. L. Root for the stationary case in “Random Signals and Noise", p, 116, McGraw 


Hill (1958). 


P(N, ty, to) = H(N, ty, ty) exp [ - h (x) dx] 
(10) 
1) A M 
H(N, ty, ty) = h (A) n@~ | h (N) dA dN 
ty 


The expression for H can be further simpified (Appendix I) and we have 


2 N 
(11) H(N, ty ty) = aT [f h (x) dx] 
ty 


By equations (10) and (11) the probability P(N, t,, to) can be written 


tg 
N 
(12) P(N, to) = | f exp]- 
ty 


It can also be verified directly by substitution that the probability (12) satisfies the 
differential equation (8) of P(N, t,, ty). 
A comparison of (12) with the Poisson distribution 


N 
(13) P(N) = efit; =meanof N 


reveals that (12) is a Poisson distribution with variable parameter 


tg 
i = ] S | h (x) dx = the mean value of h(t) in the interval (ty» ty) 


(14) 4 


| 
. 
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THE AUTOCORRELATION FUNCTION 


The autocorrelation function R(t) of a function f(t) it defined [1] as 


+T 
(15) R(t) = f(t f(t-t) dt 
-T 


In the following we shall denote by R! (7) the function which has a finite range of 
integration (-T, +T) before the limit is taken in (15), 

The current i (t) of the diode or the ionization chamber can be considered the sum 
of short pulses of charge q occurring at times t,. Assume, at first, that the pulses are 
of infinitesimal duration so that they may be represented by delta functions 4(t-t,). 
Thus, we may write 


(16) ig (t) 


By the definition, the finite autocorrelation function R’, (T) of i 6 (t) is 


11 


+T 
-T 
+T 
(17) = d(t-t) 7) dt 
om | | 
k 


The restrictions on t, and t; are due to the range of integration (-T, +T). The deltas 
on the first and second line of (17) are shifted to t, and t; + t respectively, and those 
deltas that are outside this range do not belong to the integral. 

The double sum in (17) can be divided into two parts by treating separately the 
cases k=j and k # j. 


2 2 


Because of the restriction on t, in (17), the summation in the first sum of (18) is 
over the pulses in the interval (-T,+T) only. The restrictions on t, and tj in the 
second sum are the same as in (17). As To, the first term of (18) becomes 
obviously 


2 
(19) lim & = 


where ft has the same meaning as in (14). 

When the length of 2T is increased, the number of pulses in the double sum of 
(18) grows. Finally the deltas §(1+ tj - th) become infinitely dense. Then the 
double sum can be replaced by the probability density PKj (8) of finding two pulses 
an interval © = t, - tj apart from each other. The above event occurs when the 
following two events happen together: 


a) There is a pulse between t=t, and t =t, + d@ 
b) There is a pulse between t = tj and t = tj + 


12 


These two events are statistically independent, and by equation (1) their probabilities 
are 
b) h (tj) d@ = h (t, - 6) dé = h (ty - ® dt, 


The joint probability of events a) and b) is the product of their probabilities, or 


According to equation (17), t, can be any instant in the interval (-T, +T), and the 
total probability is obtained by integrating (20) with respect to t, over the same 
range. 


+T 


(21) Pyj (0) = J h (ty) h (t, - 6) dty 


By this equation, the second part of (18) can now be written 


+T 


2 
-T 


where ae denotes the finite autocorrelation function of h (t), Returning back to the 
variable t and taking the limit T — oo of (22), the autocorrelation function of i 5 (t) 
can be written 


= + R, (7) 


(23) +T 


R, (7) = h(t)h(t- 7) dt 


-T 


Let us assume that all the current pulses are of a finite non-zero duration t, and 
otherwise identical, except that they start at different times t,. Denote the pulse 
starting at t=0 by i, (t) so that 


i 
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#0 when OStsu, 


(24) i, 
1= at other times 


The total current i (t) is the sum of the terms ip (t - t,), where i, (t) has been shifted 
to start at t,. The current i (t) can be conveniently given in any one of the following 


forms: 


+@ 
i (t) = i, (t- = ig (x) 6 (t- t, x) dx 
+@ 
= - tp dy 
(25) -@ 


= 2 b(t - 


= Gio 


where denotes the convolution *) and ig (t) has the same meaning as in (16). 
It follows from equation (25) (see Appendix II) that the autocorrelation function R, ( tT) 
of i(t) is 


+0 +0 
J ip -y)R ~y) dxdy 


-@ -@ 


R(T) 


(26) 


q 


2 
The convolution i, (tT) Wi, (- t) is often denoted by i,(7) and called the 
convolution square of i,. Using this notation and (23) we have 


+@ +@ 


x) = f f(x)g(t-x)dx= f (t - x) g(x) dx 


-@ 


4 
. 
| : 
= 
{ 
| 
i 


Ri (7) 


2 2 
(T) + (7) WR, (7) 


2 
ae (7) 


which is a general expression for the autocorrelation function of i (t) in (24). 
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THE POWER SPECTRUM 


The statistical power spectrum S (f) of the noise (f = frequency) is the average 
spectrum of the ensemble of sample functions that can be generated under the same 
conditions, The power spectrum of a particular sample function is not exactly S (f), 
but the average of a great number of sample functions approaches S (f). 

The power spectrum S (f) is known [1] to be the Fourier transform of the auto- 
correlation function R (TT). 


@ 
(28) (f) = | R(T) j=V-1 
-@ 


The Fourier transform of the convolution of two or more functions is the product 
of the transforms of each function. Thus, by equations (27) and (28) applied to R; 
and R,. 


where I, (f) is the transform of i, (t) and S; (f) and S} (f) are the power spectra of 
i(t) and h(t). 

Equation (29) is a general expression for the spectrum of a nonstationary shot 
current. In the following we shall discuss the special case of h (t) being a periodic 
function of time with a fundamental frequency fy. Then h(t) cam be given as a 
Fourier series 


3 


+0 


(30) h(t = > Cm exp (2 2 jmf,t) 
m=-@ 


If h(t) is truly periodic from - 00 to +00, the autocorrelation function is also periodic 
with the same fundamental period and [1] 


+@ 


(31) = > 


m=-@ 


2 
exp (2 %jmf, T) 


The Fourier transform of the periodic function (31) consists of delta peaks at f = mfp 
having weights |<m| 2 Hence 


+0 


(82) = > 


m=: @ 


Cm $f - 


Keeping in mind that n = C,, substitution of (32) into (29) yields finally 


2 
+ 


2 
(33) = |72 | (f - mfo) 


The first term of (33) gives the continuous noise spectrum, upon which are super- 
imposed the discrete “signal” frequencies mf, of the second term (see Fig. 1.). 
At frequencies f < 1/t,, where ty is the width of the current pulse i, in (24), 


the transform 


tw 
(34) f iy () dt=q=me ; f <ift, 
0 


where m is the number of electrons in one pulse and e is the charge of the electron. 
Thus the continuous noise spectrum has in the low frequency domain the constant 
value 


17 


y | 
B 
A 
ty 
ty ty x 
Fig. 1. 
2 2 
(35) Si = fi (me) 


The average current is obviously i =iime and the noise spectrum has the value mei, 
In the ionization chamber m is usually several hundreds, the number of ion pairs 
produced by a quantum, In the diode current the pulses consist of single electrons, 
m =1, and the noise is ei, which is known as the Schottky formula 

In the stationary case h (t) = = constant 


ii when m=0 
(36) 


0 when m #0 


and (33) becomes using (34) and putting m=1 


Ip + 6H 


(37) =i 


an expression derived earlier [1] for the stationary shot noise of a temperature- 
limited diode, 


q 
i 
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APPENDIX I 


In order to arrive at the expression (11) for H (N, ty, ty) from that given in (10) 
let us discuss separately the cases N=1, 2, 3, -. 
When N =1, by the definition (10) 


t2 
(38) H (1, th ty) = | h (x) dx 
When N =2 
to x 
(38) H(2, ty, dxdy 
We may write this double integral as 
sysx 
(40) H(2, ty, = h(x) h(y) dA; A: 
A 


where A is the triangle, shown in Fig. 2., in the xy-plane and defined by the ine- 
qualities (40). 
The function h (x) h (y) under the intergral (40) is symmetrical with respect to x 


‘ 
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Fig. 2. 


and y. Hence, the same integral over the triangle B (Fig.2.) is also equal to 
H (2, iso to), or 


tg to 
A+B hoy 
(41) 
tg 2 
=} ax 
q 


When N = 3, we have a pyramid in a threedimensional space instead of the triangle 
A of Fig. 2, Six equal pyramids, which are in symmetrical positions with respect 
to the three coordinate axes, make up a cube with an edge of length tg - t, 
and 


i 
0 1 2 3 3 5 f/f. 
q 
| 
‘ 


3 
(42) H(3, to) = f h (x) dx 
4 


By similar arguments we can conclude generally 


to N 
(43) H(N, ty, ty) = tr h (x) dx 


| 
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APPENDIX II 


By the definition (15) the autocorrelation function R, of i (f) of equation (25) 


+T 
-T 
+T 


(44) * f qio q to (t T)Wig dt 
T 


+T +0 
lim 1 


-@ 


By changing the order of the integration process and the process of taking limits we 
have heutristically 


-@ -@ 


+O +0 +T 
1 
(45) R(T) = f J is (t-x)ig (t- - y) dt} dxdy 
-T 


i 
| 
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A comparison of the expression in the brackets with the definition (15) shows that 
the former is equal to Ry (t - x +y)/q%, and 


+O +0 


When the substitution z = x - y is made and the integration over x is carried out 
keeping y constant (dx =dz), we have 


+0 +0 
1 
f f ig (X) iy (x - 2) dx] dz 
-@ -@ 


(47) 


1, 
ig (T) Wig (- (7) 


| 
| 
| 
| 
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